Fractal Descriptors in the Fourier Domain Applied to Color Texture Analysis 
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The present work proposes the development of a novel method to provide descriptors 
for colored texture images. The method consists in two steps. In the first, we apply 
a linear transform in the color space of the image aiming at highlighting spatial 
structuring relations among the color of pixels. In a second moment, we apply a 
multiscale approach to the calculus of fractal dimension based on Fourier transform. 
From this multiscale operation, we extract the descriptors used to discriminate the 
texture represented in digital images. The accuracy of the method is verified in the 
classification of two color texture datasets, by comparing the performance of the 
proposed technique to other classical and state-of-the-art methods for color texture 
analysis. The results showed an advantage of almost 3% of the proposed technique 
over the second best approach. 
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I. INTRODUCTION 



Texture is a visual attribute which is present in most of the nature images. The texture 
quantification, identification and classification is a very important problem that is defined in 
computer vision as the study of pixel patterns in an image region. Although this attribute 
is naturally processed by natural vision and easily comprehended by humans, there is no 
formal definition for it. Indeed, textures are complex visual patterns formed by arrangements 
of pixels, regions or even set of patterns formed by other visual attributes, such as shape 
or color. These patterns can be composed by completely distinct factors, such as pixel 
organization or even its disorganization. In fact, depending of the context, even the noise 
can be considered as a sort of texture. These characteristics of the texture attribute make 
it special and hard to be well defined. 

Along the last years, a lot of methods has been developed for the analysis of textures. 
Such interest in texture analysis methods may be comprehended by the richness of the 
texture attribute in images analyzed in pattern recognition problems. 

Basically, the texture analysis methods can be divided into 4 categories 4 , that is, the 
structural methods, in which the texture is described as a set of primitives well defined; 
statistical methods, in which the texture is represented through non-deterministic measures 
of distribution; spectral methods, based on the analysis in the frequency domain and model- 
based methods, based on the mathematical and physical modelling of the texture image. 

Among the model-based methods, the fractal model has presented a large projection 
recently in the description of textures in a wide number of problems 5-7 . Most of these 
methods employs the fractal dimension direct or indirectly for the representation of the 
texture. In recent years, however, a family of methods 3 ' 8-10 was developed through extracting 
a set of features from the fractal modelling, unlike the fractal dimension which is only a 
unique number. Generically speaking, these methods provide the called fractal descriptors, 
capable of representing a texture with a higher degree of richness than the simple fractal 
dimension. 

The ability of fractal features in the description of textures is related to the nature of 
fractality concept. The fractal dimension measures the complexity of a structure, which 
in its turn corresponds to physical important properties of a material, like the roughness, 
the reflectance, among others. Finally, such properties are strong stimuli in our visual 
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identification of texture images, allowing the classfication of objects based on their texture 
aspects. In this way, fractal theory becomes a worthy tool in the automation of this process. 

Actually, the use of term "fractal" in such kind of application may bring some controversy, 
once the analyzed images are not real fractals. Even because fractals are only mathemat- 
ical entities lacking any perfect representation in the real world. The interested reader 
may appreciate this discussion for example in a letter exchange involving Avnir et al. and 
Mandelbrot 11 ' 12 . Inasmuch as this debate is not in any way finished, the most acceptable 
approach in the literature is that presented in Carlin 13 . There, objects from real world are 
measured through fractal metrics even when they get away from fractal concept. In this 
case, the fractal measures acts as a complexity metric of the real world object. 

This work proposes a novel technique for the extraction of fractal descriptors based on 
the Fourier fractal dimension method 14 from colored textures. The color is an important 
attribute in texture images, mainly those extracted from natural scenes 15 . The method 
consists in representing the color image in a new space color through the linear transform 
described in Geusebroek et al. 16 . This transform aims at emphasizing the relation between 
the pixels color and their spatial distribution. Posteriorly, we apply the method for the 
calculus of fractal dimension by the Fourier transform. In this method, the dimension 
is obtained from a curve relating the power spectrum of the Fourier transform and the 
frequency. Instead of simply using the dimension value, this work proposes the use of the 
whole curve to provide the descriptors of the texture. 

The proposed method takes some important advantages over classical fractal signature 
techniques like that based on wavelets 17 or multifractal 18 . One of such advantages is that 
the technique here presented gathers information from frequency domain inherently, allowing 
the capturing of details and patterns which escapes from the conventional spatial analysis. 
Besides, the space used allows the expression of colors and spatial distribution of pixels as be- 
ing a related entity. This relation is important in many applications involving color texture. 
Moreover and not less relevant, the method shows a simple computational implementation 
and presents a low computational effort. 

The performance of the proposed technique is tested in a comparison with other classi- 
cal and state-of-the-art methods for color texture analysis, namely, chromaticity moments, 
histogram ratio and multispectral Gabor. The comparison is achieved in the classification 
of two color texture datasets. 
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This work is divided into 8 sections. The following addresses the definitions and theoret- 
ical aspects of fractal theory. The third section describes the Fourier fractal dimension. The 
fourth presents the concept of fractal descriptors. The fifth describes the proposed method. 
The following shows the experiments employed. The seventh section discusses the results 
and the last one does the conclusions. 

II. FRACTAL 

The existence of strange objects which do not obey the rules of the traditional euclid- 
ian geometry is known from mathematicians from some centuries ago. Nevertheless, the 
formalization and denomination of such objects is due to Benoit Mandebrot in the 1970 
decade 1 . 

A fractal is defined as a set whose Hausdorff-Besicovitch dimension exceeds strictly the 
topological (Euclidian) dimension. Such fact has as a consequence the fact that the fractals 
are dynamical systems with infinite complexity. Besides, fractals are self-similar, that is, 
each part of the object is a similar copy of the whole. It is noticeable that this repetition 
of patterns along different observation scales is also present in objects found in the nature, 
like the embranchment of a river, a tree or of the lung alveoli, or still in the nervures of a 
plant leaf, in a cloud, in a coastline and in many other cases 1 . 

In computational vision problems, we are interested in finding descriptors which charac- 
terize the objects in analysis. From the similarities between the aspect of natural objects 
and the objects studied in fractal theory, researchers started to study the application of a 
fractal descriptor to objects from the real world 13 . The most relevant and safe descriptor 
for this purpose is the fractal dimension. 

A. Fractal Dimension 

The literature presents several definitions of dimensions which are generally named as 
fractal dimension 19 . Among these, we can cite the Hausdorff-Besicovitch dimension, the 
packing dimension, the Renyi dimension, the box-counting dimension, etc. 

A common point among all these methods is that they are based on the idea of measuring 
at 5 scale. The analyzed object is measured for different values of 5 and at each different 
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value, the details smaller than 5 are neglected. The fractal dimension thus must express the 
behavior of the measure as 5 — > 0. In a fractal object, a measure M S (F) of a set F must 
generally obey a power law: 

M S (F) ~ c5~\ (1) 

where c is a constant and s is the fractal dimension of F. The value of s can therefore be 
obtained from: 

s = _, iffi M (2) 

5->o /o#(<5) 

In a general way, Mg(F) must be a homogeneous function with degree d yielding to the 
power law 

M S (F) ~ c5 d " s , (3) 

Here, we describe briefly the development of the first and perhaps most important measure 
of the fractal dimension, e.g., the Hausdorff dimension. For a general set F G 9ft n , the 
Hausdorff dimension is defined by the following expression: 

dim H (F) = {s}\ inf {s : H S (F) = 0} = sup {H S (F) = oo} , (4) 

where H S (F) is the s-dimensional measure of F, defined through: 

H s (F)=]imH s e (F), (5) 

where 

H s e (F) = inf |^ \Ui\ s : U t is an e-cover of f| . (6) 
III. FOURIER FRACTAL DIMENSION 



The literature still shows alternative definitions for the fractal dimension. Among these 
definitions, one of the most important is the Fourier fractal dimension 19 . For the calculus of 
this dimension we must initially define the Fourier transform of a mass distribution /i e 5J n , 
that is, a measure on a bounded subset of 9ft n such that /i(5i n ) is positive and finite. The 
transforms are defined by: 

= f e ix ' u dfi(x), (7) 
where u is a generic subset of !R n and a; is the Fourier space counterpart of u. 
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In the following we employ an analogy from the classical mechanics for the definition of 
the fractal dimension. We use the concept of s-potential of a mass distribution \i over a 
point x in 9ft n , given by: 

Ps{x) = I 1 , s d/i(y), (8) 
J \x-y\ 8 

where y is an auxiliary variable. 

By still extending the Physics analogy, the potential energy c s may be obtained through: 

e s (/i) = {2n) n c J Z(p s )(u)Z(tI(uj)du, (9) 

in which c is a constant dependent on s and n and x is the complex conjugate of x. In this 
way: 

e s (/i) = (2vr) n c J \u\ s - n \1(fi(u))\ 2 du. (10) 

From a theorem developed in Falconer 19 , if there is a mass distribution fi(u) on the set 
S G 3? ra for which the expression 10 is finite for some value(s) of s, so the Hausdorff dimension 
of S has its lower limit in s. Particularly, if \%(fi(u))\ < 6|w|~*/ 2 , for a constant value b, then 
e s (/x) always converges if s < t. The greatest t for which there is a mass distribution /i on S 
is called the Fourier fractal dimension of S. 

For practical purposes, the method for the calculus of the Fourier fractal dimension 
described in Russ 14 is applied to a gray-scale (intensity) image /. 

In this case, we have a 2D real- valued image with size N x N and the Fourier 

transform / is expressed through: 

N-l N-l 

I(u,v) = 5>(M)exp-^ + ^, (11) 

1=0 jr=0 

where u and v are respectively the horizontal and vertical frequency. The total frequency 
/ is given by / = y/u 2 + v 2 . Another important measure is the power spectrum P, given 
through P = P. 

Russ 14 demonstrates that there is an exponential relation between the frequency / and 
the power spectrum P in the Fourier spectrum of /: 

p oc r a , (12) 



He still affirms that the exponentiation parameter a may be used for the estimation of the 
fractal dimension D of the texture in /. In the practice, a is calculated as being the slope 
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of the curve log(P) x log(f). The dimension is easily estimated through: 



In digital images applications, the Fourier transform is calculated by classical optimized 
techniques, like Fast Fourier Trasnform 20 . The Fourier spectrum is divided into radial rings 
(corresponding to frequency bands). Thus, the variable P in the previous expression corre- 
sponds to the power spectrum averaged over each ring and the frequency / corresponds to 
the average distance of each ring from the center of the spectrum. The Figure 1 illustrates 
the process. The Figure 2 illustrates the dimension calculated by this practical method. 




Figure 1. Fourier spectrum is divided into radial rings and the power spectrum P is averaged over 
each ring. The frequency / corresponds to the average distance of each ring to the center of the 
spectrum. 



IV. FRACTAL DESCRIPTORS 

Although the fractal dimension is a good descriptor for textures, shapes, contours, etc., 
they become inefficient in tasks which require a greater precision in the description of the 
object. Fractal descriptors have arrisen with the aim of filling this gap and provide a more 
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Figure 2. Calculus of the fractal dimension of textures, (a) Original texture, (b) Power spectrum, 
(c) Log-log curver of power spectrum x frequency. 

precise technique for the characterization of the image. The Figure 3 shows graphically the 
importance of fractal descriptors. 

The first known work applying the concept of fractal descriptor is Manoel et al. 2 which 
uses the technique named Multiscale Fractal Dimension (MFD) to obtain the descriptors. 
In this approach, the fractal dimension from the object is inferred at different observation 
scales. In Manoel et al. 2 , the authors used the Minkowski sausage method for the calculus 
of the fractal dimension. In this method, the object is dilated by a variable radius r and 
the area of the dilated object (number of pixels) is called the dilation area A{r). The 
fractal dimension is estimated from the slope of the curve log(A(r))) x log(r). Instead of 
simply obtaining the dimension, the authors used the whole log(A(r)) curve to compose the 
descriptors by measures extracted from the curve, like the peaks, the area under the curve, 
etc. 

In its turn, in Plotze at al. 9 and Bruno et al. 3 we have an application of MFD in which the 
derivative of log(A(r)) is used in order to provide the fractal descriptors. Bruno et al. 8 still 
applies the MFD to the analysis of textures, mapped onto surfaces and using the volumetric 
Minkowski ge method. 

V. PROPOSED METHOD 

This work proposes a novel method for the extraction of fractal descriptors from colored 
texture images. 

A lot of methods has been described in the literature for the extraction of features from 
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Figure 3. Illustration of the richness of texture fractal descriptors. At left, two textures with 
similar fractal dimensions. At right, the fractal descriptors for each texture and the clear visual 
distinction between them. 



colored textures with the aim of solving problems like classification and segmentation in 
different application fields 21-23 . However, many of such works do not take into account the 
spatiality of the color, that is, the relation between the color of a pixel and its position in 
the image or in a specific neghborhood. 

In order to deal with this situation, Geusebroek et al. 16 proposed an interesting method 
based on physical characteristics of colors. Roughly speaking, the method consists in a linear 
transform from the original color space into another physical space. In the special case in 
which the original space is RGB (Red-Green-Blue) like it is our case, the transformation is 
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represented by the simple expression: 



i Ex ^ 

y E\\\ j 



( 



\ 



0.06 0.63 0.31 
0.19 0.18 -0.37 
0.22 -0.44 0.06 



G 



(14) 



where R, G and B are the original color channels and E\, E X \ and E X \\ are the transformed 
channels. E\ corresponds to the convolution of color energy (wavelength) with a gaussian. 
E\\ represents the same convolution with the first derivative of a gaussian, while E\\\ is 
the same with the second derivative of gaussian filter. In Hoang et al. 24 , the authors apply 
the classical Gabor filters to the transformed color space obtaining interesting results. 

In this work, initially we apply the transform described to the texture image. In the 
following, for each one of the channels E X) E X \ and E\ X \ we extract fractal descriptors 
based on the Fourier fractal dimension described in the Section III. For the calculus of these 
descriptors, instead of simply calculate the fractal dimension by the Equation 13, we use all 
the values of log(P) in the curve, that is, the logarithm of the whole power spectrum. Thus, 
we appy a multiscale transform to the curve, with the aim of capturing the fractal behavior 
at different observation scales, in a similar manner to that described in Plotze et al. 9 . 

Essentially, a multiscale transform is a mapping from the original signal u(t) onto a 
function U (b, a), where b is related to the original variable t and a is the scale parameter. The 
literature presents several approaches for the calculus of multiscale transform 25 ' 26 . Based on 
empirical results we opted for the use of space-scale approach. In this solution, the transform 
is formally represented through: 



{(b,a)\a,be 5ft, a > 0,6 e {U\t, a)} zc }, 



(15) 



in which . zc represents the zero-crossings of . and U 1 (t,a) expresses the convolution of u(t) 
with the first derivative of the gaussian g\, given by: 



U 1 (t,a)=u(t)*g 1 a (t), 



(16) 



where a denotes the smoothing parameter of gaussian, referenced in the most of textbooks as 
a. The best results were achieved by calculating the derivative through the Fourier property 
and projecting U(b, a) onto the axis corresponding to a = 0. This multiscale representation 
based on Equations 15 and 16 uses the gaussian kernel concept. It states, using a result 
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from partial derivative theory, that the gaussian filter associated to the first derivative is the 
unique operation capable of representing an image or signal under different scales without 
adding any spurious element in the process. More details may be seen in Witkin 25 

Finally, we concatenate the fractal descriptors from each channel, generating the final 
color Fourier fractal descriptor. The steps of the method are depicted in the Figure 4, while 
the Figure 5 summarizes the algorithm process. 

By joining the power of fractal theory and more specifically fractal descriptors in the 
description of natural textures and the efficiency of the color approach described in Hoang 
et al. 24 we obtain a powerful descriptor for natural colored textures. Such descriptor is 
capable of capturing complex patterns in the texture which are capital for a complete and 
precise identification of a real scene. The Figure 6 shows in a simple example the potentiality 
of texture discrimination present in the proposed descriptor which may be seen even visually. 



VI. EXPERIMENTS 

The performance of the proposed technique is tested by the classification of samples 
from VisTex 27 , a classical dataset of colored textures and USPTex, a dataset developed 
in the research group of the authors, which is composed by images of natural textures, 
photographed in high resolution. The Figures 7 and 8 show some image samples from each 
dataset. 

The dimensionality of proposed descriptors is proportional to the dimension of the image. 
Thus, for the VisTex dataset, the signature are composed by 120 descriptors, while in 
USPTex we use 132 descriptors. We have not used other resolutions for descriptors, once 
such approach may compromise the fractality measure, in which all multiscale levels has 
equal importance for the description of the texture. 

The results obtained were compared to the use of other classical methods for descriptors of 
colored textures, that is, color Gabor 24 , color histogram ratio 28 and chromaticity moments 22 . 
The descriptors are classified by the well known K-Nearest Neighbor (KNN) method 29 , with 
k — 1 (empirically determined) and using a 10-fold cross-validation process. The comparison 
is done in terms of the correctness rate with its confidence interval and the confusion matrices 
for each descriptor and dataset. 
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Figure 4. A scheme of the proposed method. From up to down, the original texture, the transformed 
channels, the Fourier curve for each channel and the final descriptor. 
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1 Load image 

2 Extract E_lambda channells using the transform matrix 

3 For each channell 

4 Calculate Fourier tranform 

5 Divide spectrum into radial rings (frequency band) 

6 For each ring 

7 power spectrum = squared magnitude 

8 end 

9 Fractal descriptors = multiscale (log (power spectrum ) xlog ( frequency ) ) 

10 end 

11 Final descriptors = concatenate ( fractal descriptors) 



Figure 5. Proposed method generic algorithm. 
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Figure 6. The ability of the proposed descriptors in the discrimination of texture classes. Left, 
we see texture images from two classes and their respective descriptors. Right, the descriptors are 
plotted in a same graph, showing visually the high discrimination potential. 
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Figure 7. One image sample from each class of the Vistex dataset. From up to down, left to right: 
Sand, Tile, Water, Bark, Fabric, Food and Metal. 

VII. RESULTS 

In a first moment, we tested the classification performance of some metrics extracted from 
MFD curve, as it is done in the original MFD work 2 . Here, we employed some fractal and 
statistical metrics, e.g., fractal dimension of channells 1, 2 and 3, mean, standard deviation, 
kurtosis, skewness, second and third order moments and the combination of these measures. 
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Figure 8. Some image samples (one from each class) from the USPTex dataset. 

The Table I shows the correctness rate and associated error for each approach in Vistex 
dataset. Thus, we see that the use of whole MFD curve in the composition of descriptors 
provided the best result. From now, we show results for the use of the whole curve in the 
tested datasets. 

Initially, we see in the Table II the global correctness rate for each compared descriptor 
for Vistex and USPTex dataset. It is clearly noticeable that the proposed Fourier method 
presented the best result in the classification of both datasets. The proposed technique 
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Table I. Percentage correctness rate and respective confidence interval in the classification of Vistex 
dataset using statistical metrics extracted from MFD curve. 

Metric Correctness rate 



FD Channell 1 30. 


.51 


± 


0.20 


FD Channell 2 33 


.11 


± 


0.23 


FD Channell 3 35. 


.38 


± 


0.14 


Mean 


32. 


.46 


± 


0.21 


Std. Dev. 


26. 


.62 


± 


0.17 


Kurtosis 


20. 


.12 


± 


0.27 


Skewness 


29. 


.54 


± 


0.39 


2 th Moment 


26. 


.62 


± 


0.30 


3 th Moment 


29. 


.87 


± 


0.40 


Combined 


93. 


83 


± 


0.21 


Whole Curve 


95, 


12 


± 


0.10 



presented an advantage of 2.8% in Vistex dataset and 2.7% in USPTex over Gabor method, 
the second best technique in this experiment. 

Table II. Percentage correctness rate and respective confidence interval in the classification of the 
tested datasets by the compared descriptors. 

Method Vistex USPTex 

Moment 68.83 ± 0.33 32.06 ± 0.05 
Histogram 78.89 ± 0.21 41.49 ± 0.18 
Gabor 92.53 ± 0.18 86.47 ± 0.04 
Fourier 95.12 ± 0.10 88.83 ± 0.07 



In the Table III we present the confusion matrices for each descriptor method. In this 
matrix, each raw (or column) represents a class and the value in raw i and column j expresses 
the number of objects of class i, but classified as being from class j. The ideal method (with 
a 100% correctness rate) must present a diagonal confusion matrix. 

For the Vistex dataset, the numbers correspond to the following classes (exemplified in 
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the Figure 7): 1-Bark, 2-Fabric, 3-Food, 4-Metal, 5-Sand, 6-Tile and 7- Water. As expected, 
the greater the correctness rate, more the confusion matrix presents diagonal aspect. Par- 
ticularly, the proposed descriptor presented its best performance in the classes 3, 4 and 
5. Good results are also observed in classes 1 and 6. Specially, in class 6, the proposed 
technique presented a relevant advantage over Gabor method. Gabor misclassified elements 
from class 6 as being from class 1, 2 and 5. This is explained by the self-similarity present 
in these classes. The fractal method, as expected, captured more faithfully the self-similar 
nuances. The confusion matrices from USPTex were represented in a different graphical 
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(c) (d) 
Table III. Confusion matrices for the classification of Vistex dataset using the compared descriptors, 
(a) Chromaticity moment, (b) Histogram, (c) Gabor. (d) Fourier. 

manner, by using surface figures in the Figure 9. In this representation, each position in the 
matrix is represented by a surface point and the height of the point, according to the legend 
on the axis, determines the value in that position. Looking at each figure, we observe that 
the matrix from the Fourier method presented a more continuous diagonal with the highest 
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points, justifying the greater number of samples correctly classified. The Gabor method 
presented bad results around the class 110. 




Expected class Expected class 

(c) (d) 
Figure 9. Surface visualization of confusion matrices for the classification of USPTex dataset using 
the compared descriptors, (a) Chromaticity moment, (b) Histogram, (c) Gabor. (d) Fourier. 



VIII. CONCLUSION 

This work proposed a novel technique for the calculus of descriptors for colored textures. 
The method uses the Fourier spectral dimension associated to the spatial color transform 
proposed in Geusebroek et al. 16 . Initially, the transform is applied to the original texture 
images in RGB space. Following, the classical Fourier transform is applied to the image and 
the values in the curve log (power spectrum) x log( frequency) are used as descriptors for the 
texture. 

The accuracy of the proposed method was verified by a comparison with other classical 
techniques in color texture analysis. The results demonstrated the power of the novel tech- 
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nique once Fourier descriptors presented the greater precision in the classification of two 
complex datasets. Specially, the proposed technique showed great efficiency in capturing 
self-similar patterns in the texture. 

The results suggest strongly that color Fourier fractal descriptors are an interesting al- 
ternative to be used in the solution of problems in which the description of an object is 
primordial, like tasks involving segmentation and classification of objects represented by 
their texture. 

IX. ACKNOWLEDGEMENTS 

Odemir M. Bruno gratefully acknowledges the financial support of CNPq (National 
Council for Scientific and Technological Development, Brazil) (Grant #308449/2010-0 and 
#473893/2010-0) and FAPESP (The State of Sao Paulo Research Foundation) (Grant # 
2011/01523-1). Joao B. Florindo is grateful to CNPq(National Council for Scientific and 
Technological Development, Brazil) for his doctorate grant. 



19 



REFERENCES 



1 B. B. Mandelbrot, The Fractal Geometry of Nature (Freeman, NY, USA, 1975). 

2 E. T. M. Manoel, L. da Fontoura Costa, J. Streicher, and G. B. Miiller, "Multiscale 

fractal characterization of three-dimensional gene expression data," in SIBGRAPI (IEEE 

Computer Society, 2002) pp. 269-274. 

3 0. M. Bruno, R. de Oliveira Plotze, M. Falvo, and M. de Castro, "Fractal dimension 
applied to plant identification," Information Sciences 178, 2722-2733 (2008). 

4 A. Materka, M. Strzelecki, T. Analysis, M. A. Review, A. Materka, and M. Strzelecki, 
"Texture analysis methods - a review," Tech. Rep. (Institute of Electronics, Technical 
University of Lodz, 1998). 

5 R. Quevedo, F. Mendoza, J. M. Aguilera, J. Chanona, and G. Gutierrez-Lopez, "Determi- 
nation of Senescent Spotting in Banana (Musa cavendish) Using Fractal Texture Fourier 
Image," Journal of Food Engineering 84, 509-515 (2008). 

6 L. Tian-Gang, S. Wang, and N. Zhao, "Fractal Research of Pathological Tissue Images," 
Computerized Medical Imaging and Graphics 31, 665-671 (2007). 

7 H. Millan and M. Gonzalez-Posada, "Modelling Soil Water Retention Scaling. Comparison 
of a Classical Fractal Model with a Piecewise Approach," GEODERMA 125, 25-38 (2005). 

8 A. R. Backes, D. Casanova, and O. M. Bruno, "Plant leaf identification based on vol- 
umetric fractal dimension," International Journal of Pattern Recognition and Artificial 
Intelligence (IJPRAI) 23, 1145-1160 (2009). 

9 R. O. Plotze, J. G. Padua, M. Falvo, M. L. C. Vieira, G. C. X. Oliveira, and O. M. Bruno, 
"Leaf shape analysis by the multiscale minkowski fractal dimension, a new morphometric 
method: a study in passiflora 1. (passifloraceae)," Canadian Journal of Botany- Revue 
Canadienne de Botanique 83, 287-301 (2005). 

°J. B. Florindo, M. De Castro, and O. M. Bruno, "Enhancing Multiscale Fractal Descriptors 
Using Functional Data Analysis," International Journal of Bifurcation and Chaos 20, 3443- 
3460 (2010). 

1 D. Avnir, O. Biham, D. Lidar, and O. Malcai, "Applied Mathematics: Is the Geometry 
of Nature Fractal?" Science 279, 39-40 (1998). 

2 B. B. Mandelbrot, "Is nature fractal?" Science. 279, 783-785 (1998). 



20 



13 M. Carlin, "Measuring the complexity of non-fractal shapes by a fractal method," Pattern 

Recognition Letters 21, 1013-1017 (2000). 
14 J. C. Russ, Fractal Surfaces (Plenum Press, New York, 1994). 

15 B. Manjunath, J. Ohm, V. Vasudevan, and A. Yamada, "Color and Texture Descriptors," 

IEEE Transactions on Circuits and Systems for video Technology 11, 703-715 (2001). 
16 J.-M. Geusebroek, R. van den Boomgaard, A. W. M. Smeulders, and A. Dev, "Color and 

scale: The spatial structure of color images," in Sixth Europian Conference on Computer 

Vision (ECCV (Springer, 2000) pp. 331-341. 
17 K. Muneeswaran, L. Ganesan, S. Arumugam, and K. R. Soundar, "Texture classification 

with combined rotation and scale invariant wavelet features," Pattern Recogn. 38, 1495- 

1506 (2005). 

18 D. Harte, Multifractals: theory and applications (Chapman and Hall/CRC, 2001). 
19 K. J. Falconer, The Geometry of Fractal Sets (Cambridge University Press, New York, 
NY, USA, 1986). 

20 O. Brigham, The Fast Fourier Transform (Prentice-Hall, NY, USA, 1974). 

21 T. Caelli and D. Reye, "On the classification of image regions by colour, texture and 
shape," Pattern Recognition 26, 461 - 470 (1993). 

22 G. Paschos, "Chromatic correlation features for texture recognition," Pattern Recognition 
Letters 19, 643 - 650 (1998). 

23 A. Jain and G. Healey, "A Multiscale Representation Including Opponent Color Features 
for Texture Recognition," IEEE Transactions on Image Processing 7, 124-128 (1998). 

24 M. A. Hoang, J.-M. Geusebroek, and A. W. Smeulders, "Color texture measurement and 
segmentation," Signal Processing 85, 265 - 275 (2005). 

25 A. P. Witkin, "Scale space filtering: a new approach to multi-scale descriptions," in Pro- 
ceedings..., ICASSP - IEEE International Conference on Acoustics, Speech, and Signal 
Processing (GRETSI, Saint Martin d'Hres, France, 2003) pp. 79-95. 

26 L. da F. Costa and R. M. Cesar, Jr., Shape Analysis and Classification: Theory and 
Practice (CRC Press, 2000). 

27 MIT, "Mit vistex texture database," (2010). 

28 G. Paschos and M. Petrou, "Histogram ratio features for color texture classification," 
Pattern Recognition Letters 24, 309 - 314 (2003). 



21 



29 R. O. Duda and P. E. Hart, Pattern Classification and Scene Analysis (Wiley, New York, 
1973). 



22 



